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Motivated by applications in computational anatomy, we consider a second-order 
problem in the calculus of variations on object manifolds that are acted upon by Lie 
groups of smooth invertible transformations. This problem leads to solution curves 
known as Riemannian cubics on object manifolds that are endowed with normal metrics. 
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The prime examples of such object manifolds are the symmetric spaces. We characterize 
the class of cubics on object manifolds that can be lifted horizontally to cubics on the 
group of transformations. Conversely, we show that certain types of non-horizontal 
geodesies on the group of transformations project to cubics. Finally, we apply second- 
order Lagrange-Poincare reduction to the problem of Riemannian cubics on the group 
of transformations. This leads to a reduced form of the equations that reveals the 
obstruction for the projection of a cubic on a transformation group to again be a cubic 
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1 Introduction 

In this section, we summarize the main content of the paper, motivated by potential apphca- 
tions of first-order and higher-order trajectory planning problems in computational anatomy. 

1.1 General Background 

Geodesic matching in computational anatomy. The new science of computational 
anatomy (CA) is concerned with quantitative comparisons of shape, in particular the shapes 
of organs in the human body [CMOS]. In CA, shapes are defined by spatial distributions of 
various types of geometric data structures, such as points (landmarks), spatially embedded 
curves or surfaces (boundaries), intensity, or density (regions), or tensors that encode local 
orientation of muscle fibers, etc. A fruitful approach in this burgeoning field applies the 
large deformation matching (LDM) method. In the LDM method, shapes are compared by 
measuring the relative deformation required to match one shape to another under the action 
of the diffeomorphism group [DGM98, Tro98]. This approach follows D'Arcy Thompson's 
inspired proposal to study comparisons of shapes by transforming one shape into another 
[Tho42]. More specifically, it follows the Grenander deformable template paradigm [Gre93]. 

If a Lie group G is equipped with an invariant metric, then its action on a smooth manifold 
Q induces a metric on Q called the normal metric. If the diffeomorphism group is equipped 
with a right-invariant metric, then the normal metric associated with the group action is a 
natural choice for the metric on shape space. For discussions of normal metrics induced by 
actions of the diffeomorphism group on smooth manifolds in the context of computational 
anatomy, see e.g. [YAM09, YoulO]. 

The objective in computational anatomy is to quantify the distance between two given 
shapes by computing the length of the geodesic path between them, with respect to the 
normal metric on shape space. This is equivalent to computing a (horizontal) geodesic path 
on the group of diffeomorphisms that carries one shape into the other. This horizontal 
lifting property of geodesies has been crucial in the understanding and the numerics of LDM, 
and it has clarified the close connection between ideal fluid mechanics and image registration. 
Namely, geodesic flows on diffeomorphism groups are described by Euler-Poincare equations, 
called EPDiff equations in the diffeomorphism context [HMR98, MR03, YAM09]. Geodesic 
flows on the subgroup of volume-preserving diffeomorphisms were famously identified with 
Euler's equations for ideal incompressible fluid flow in [Arn66]. Paper [MTY06] studies 
momentum conservation properties of geodesic EPDiff flows in the context of image matching. 
In particular, these geodesic flows are encoded by their initial momenta, and horizontality 
means that only momenta of a specific form are permitted. For example, landmark-based 
geodesic image matching naturally summons the singular momenta that were introduced as 
solitons for shallow water waves on the real line in [CH93] and then characterized as singular 
momentum maps in any number of dimensions in [HM04]. We refer to [HRTY04, MTY06, 
CHIO, YAM09, Via09, BGBHRll, GBRll] for further details. 

In addition to their importance in computational anatomy, normal metrics and their 
geodesies also appear in a range of problems in mechanics and control theory, see for example 
[BC96, HM04, NA05]. 
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Trajectory planning and longitudinal studies in computational anatomy. Longitu- 
dinal studies in computational anatomy seek to determine a path that interpolates optimally 
through a time-ordered series of images, or shapes. This sort of task is also familiar in data 
assimilation. Depending on the specific application, the interpolant will be required to have 
a certain degree of spatiotemporal smoothness. For example, the pairwise geodesic matching 
procedure can be extended to piecewise-geodesic interpolation through several shapes, as in 
[BK08, DPT+09]. If a higher degree of smoothness is required, then one must investigate 
higher-order interpolation methods. 

Higher-order interpolation methods on finite-dimensional spaces have been studied ex- 
tensively in the context of trajectory planning in aeronautics, robotics, computer graphics, 
biomechanics and air traffic control. In particular, the study of Riemannian cubics in mani- 
folds with curvature and their higher-order generalizations originated in [GK85], [NHP89] and 
[CS95]. Riemannian cubics are solutions of Euler-Lagrange equations for a certain second- 
order variational problem in a finite-dimensional connected Riemannian manifold, to find 
a curve that interpolates between two points with given initial and final velocities, subject 
to minimal mean-square covariant acceleration. The mathematical theory of Riemannian 
cubics was subsequently developed in a series of papers including [CSC95, CCS98, CSCOl, 
Noa03b, Noa04, Noa06a, GGP02, Kra05] and Part 1 of the current study [GBHM+10]. The 
last reference treats Riemannian cubics on Lie groups by Euler-Poincare reduction. En- 
gineering applications are discussed in [PR97, ZKG98, HB04b, HB04a, Noa03a], amongst 
others. Related higher-order interpolation methods have been studied, for example, in 
[BKOO, HP04, MS04, NP05, Noa06b, MSKIO]. We refer to [Pop07, MSKIO] and [Noa06b] for 
extensive references and historical discussions concerning Riemannian cubics, their higher- 
order generalizations, and related higher-order interpolation methods. 

1.2 Motivation 

In this paper, we consider the problem of Riemannian cubics for normal metrics, focusing 
on their lifting and projection properties. Recall that in the context of normal metrics one 
considers a manifold of objects, or shapes, that are acted upon by a Lie group of transforma- 
tions. Two distinct interpolation strategies offer themselves. First, one may choose to define 
a variational principle on the Lie group, or indeed its Lie algebra, and find an optimal path 
g{t) that transforms the initial shape q as q{t) = g{t)q, such that q{t) passes through the 
prescribed configurations. This type of higher-order interpolation was proposed in regard 
to applications in Computational Anatomy, in [GBHM+10] and studied there in detail in 
the finite-dimensional setting. Alternatively, one may define a variational principle on shape 
space itself and find an optimal curve that interpolates the given shapes. This was the ap- 
proach of [VTIO], where interpolation by Riemannian cubics on shape space was proposed, 
and existence results in the case of landmarks were given. The particular cost functionals 
that interest us here are, on the group. 
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Hamilton's principle, 5S = 0, leads to cubics in the respective manifolds G and Q. One is 
therefore naturally led to the question of how Riemannian cubics on the group of transfor- 
mations are related to those on the object manifold. This is a new question in geometric 
mechanics and its answer is potentially important in applications of computational anatomy. 

In this paper, we begin the investigation of this question. We first analyze horizontal 
lifts of cubics on the object manifold to the group of transformations. In the context of 
symmetric spaces, we completely characterize the class of cubics on the object manifold that 
can be lifted horizontally to cubics on the group of transformations. For rank-one symmetric 
spaces this, selects geodesies composed with cubic polynomials in time. We then study non- 
horizontal curves in G. We show that certain types of non-horizontal geodesies project to 
cubics in Q. Finally, we present the theory of second-order Lagrange-Poincare reduction 
for Riemannian cubics in the group of transformations. The reduced form of the equations 
reveals the obstruction for such a cubic to project to a cubic on the object manifold. 

1.3 Main content of the paper 

The main content of the paper may be summarized as follows: 

In Section 2 we outline the geometric setting for the present investigation of Riemannian 
cubics for normal metrics and their relation to Riemannian cubics on the Lie group 
of transformations. In particular, we summarize the definition of higher-order tangent 
bundles by following [CMROl]. Then we define normal metrics and recall that the 
projection G — ?■ Q which maps an element of the group to the transformed image of a 
reference shape is a Riemannian submersion. 

In Section 3 we provide the key expressions for covariant derivatives of curves, and vector 
fields along curves, both in Lie groups and in object manifolds with normal metric. 
The horizontal generator of a curve in the object manifold is introduced and expressed 
in terms of the momentum map of the cotangent lifted action. 

In Section 4 we derive the equations of Riemannian cubics for normal metrics. For ease of 
exposition we first consider a more general context and then specialize to the case of 
Riemannian cubics. Here the horizontal generator plays a crucial role. Invariant metrics 
on Lie groups are a simple example of normal metrics, as are the metrics on symmetric 
spaces. These examples are worked out in detail. We also recall from [GBHM+10] 
how Riemannian cubics on Lie groups can be treated equivalently by Euler-Poincare 
reduction. Our derivation of the Euler-Lagrange equations bypasses any mention of 
curvature. Therefore, these equations can also be used to compute curvatures by means 
of the general equation for Riemannian cubics derived in [NHP89, CS95]. This is 
demonstrated by two simple examples. 
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In Section 5 we study horizontal lifting properties of Riemannian cubics. Our form of the 
Euler-Lagrange equations is particularly well suited for this task, due to the appearance 
of the horizontal generator of curves. We characterize the cubics in symmetric spaces 
that can be lifted horizontally to cubics in the group of isometrics. We then proceed 
to the more general situation of a Riemannian submersion and state necessary and 
sufficient conditions under which a cubic on the object manifold lifts horizontally to a 
cubic on the Lie group of transformations. 

In section 6 we extend the previous considerations to include non-horizontal curves on the 
Lie group. We show that certain non-horizontal geodesies on the group of transfor- 
mations project to cubics on the object manifold. We then reduce the Riemannian 
cubic variational problem on the group by the isotropy subgroup of a reference object. 
To achieve this, we use higher-order Lagrange-Poincare reduction [CMROl, GBHRll]. 
The reduced Lagrangian couples horizontal and vertical parts of the motion, and this 
explains the absence of a general horizontal lifting property for cubics. Namely, the 
reduced equations that describe Riemannian cubics on the Lie group contain the equa- 
tion that characterizes Riemannian cubics on the object manifold, plus extra terms. 
These extra terms represent the obstruction for a cubic on the Lie group to project to a 
cubic on the object manifold. In this sense, these reduced equations fully describe the 
relation between cubics on the Lie group and cubics on the object manifold. They also 
lend themselves to further study of the questions investigated in the present paper. 



2 Geometric setting 

This section introduces the geometric ingredients used in the paper. In particular, it reviews 
Hamilton's principle on higher-order tangent bundles. It also introduces Riemannian cubics 
and some of their generalizations. It then defines normal metrics and explains how the 
projection G ^ Q mapping an element of the group to the transformed image of a reference 
shape a G Q is a Riemannian submersion. 

2.1 k*'^-order tangent bundles 

The fc^'^-order tangent bundle Tq'^ : T'^^^Q — )■ Q of a smooth manifold Q is defined as a 
set of equivalence classes of curves, as follows: Two curves q2{t) are equivalent, if and 

only if their time derivatives at t = up to order k coincide in any local chart. That is, 
gf^(O) = 5'2^(0), for Z = 0, . . . , fc. The equivalence class of a given curve q{t) is denoted by 
['?]q^o)- ^'^^ then defines T^^^'Q as the set of equivalence classes of curves, with projection 

t^q'^-.T^'^^Q^Q, [q]%^q{0). (2.1) 

Finally, the inverse image of go ^ Q by Tq^ will be denoted by Tq^^'Q. This is the set of 
equivalence classes of curves based at go- Note that T^^^Q = Q and T^^^Q = TQ. 
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Given a curve q{t) one defines the fc*^-order tangent element at time t to be 

■■= , where h : r ^ qit + r) . (2.2) 

We will sometimes use the coordinate notation {q{t) , q{t) , . . . , q^^'' {t)) to denote For 
more information on higher-order tangent bundles see [CMROl]. 

A smooth map f : M ^ N induces a map between fc^'^-order tangent bundles, 

TW/ : T(^)M ^ T(^)iV , [gJ^J^ ^ [f o q]fl^^ . (2.3) 

Therefore, a group action ^ : G x Q Q on the base manifold lifts to a group action on the 
fc*^-order tangent bundle, 

$(^) : a X T(^)g ^ T(^')Q , $f : [q]% ^ [% o . (2.4) 



2.2 k^'^-order Euler-Lagrange equations 

A fc*''-order Lagrangian is a function L : T^'^^Q — )• M. In the higher-order generalization of 
Hamilton's principle, one seeks a critical point of the functional 

J[q] := f L{q{t)Ait).....q^'\t)) dt (2.5) 
Jo 

with respect to variations of the curve q{t) satisfying fixed end point conditions g'^')(0) = q^'' 
and g*-'^(l) = qf^ for / = 0, . . . , A; — 1. A curve q{t) respecting the end point conditions satisfies 
Hamilton's principle if and only if it satisfies the k^^- order Euler-Lagrange equations 

We will use from now on the standard (5-notation for variations. Let e t— )■ /e be a variation of 
a quantity f = fo- Define 

£=0 

Hamilton's principle, for example, then takes the simple form 6 J = 0. 



Examples: Riemannian cubic polynomials and generalizations. Riemannian cu- 
bics, as introduced in [viK85, NHP89] and [CS05], generalize cubic polynomials in Euclidean 
space to Riemannian manifolds. Let (Q, 7) be a Riemannian manifold and denote by the 
covariant derivative with respect to the Levi-Civita connection V for the metric 7. Let || • || 
be the norm induced by 7. Consider Hamilton's principle (2.5) for k = 2 with Lagrangian 
L : T^^^Q R given by 

(2.8) 



L{q,q,q) 



D 
Dt' 
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This Lagrangian is indeed well-defined on the second-order tangent bundle T^'^^Q, since in 
coordinates 

§-/ = q' + r^,{qWq\ (2.9) 

where T^j{q) are the Christoffel symbols at the point q. Denoting by R the curvature tensor 
defined by R{X, Y)Z := Vx^yZ - Vy^xZ - V[x,y]Z for any vector fields X,Y,Z e X{Q), 
the Euler-Lagrange equation is 



f D \ 

— + g(t),g(t)J q{t) 



0. 



(2.10) 



A solution of this equation is called a Riemannian cubic, or cubic for short. These are the 
curves we shall study in this paper. 

We also mention two generalizations of Riemannian cubics. The first one consists of the 
class of geometric k-splines [CSC95] for k > 2 with Lagrangian L : T^'^^Q — )■ R, 



L{q,q,...,q^''>) 



D 



k-l 



Dt^~^ 



(2.11) 



Note that the case k = 2 recovers the Riemannian cubics. The Euler-Lagrange equations are 
[CSC95] 

" (2.12) 



j=2 ^ ^ 

The second generalization comprises the class of cubics in tension; see, for example, [HB04b]. 



2.3 Normal metrics 

Group actions. Let G be a Lie group with Lie algebra g, acting from the left on a smooth 
manifold Q. We denote the action by 



$:Gxg^g, {g,q)^gq:=<^g{q). 



(2.13) 



The infinitesimal generator of the action corresponding to ^ G is the vector field on Q given 
by 



d 
dt 



exp{t^)q. 



(2.14) 



t=o 



In accordance with (2.4), the tangent lift of $ is defined as the action of G on TQ, 

GxTQ^ TQ, (g, v^) ^ gv, := T^g{v,), (2.15) 

with infinitesimal generator E,tq corresponding to ,^ G 0. Note that we have the relation 

TTQ{iTQ{v,)) = ^Q{q\ (2.16) 

where tq : TQ — > Q is the tangent bundle projection. Similarly, one defines the cotangent 
lifted action as 

GxT*Q^ T*Q, (g, a,) ^ ga, := (TVO*(«,)- (2.17) 
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The momentum map 3 : T*Q ^ q* associated with the cotangent hft of $ is determined by 

(J ("<?)> Og.xg = K> ^QiQ))T'QxTQ ' (2-18) 

for arbitrary aq E T*Q and ^ G 5. 

Normal metrics. Let G be a Lie group acting transitively from the left on a smooth 
manifold Q. Let 7^ be a right-invariant Riemannian metric on G. We will now use the 
action of G on Q in order to induce a metric 7q on Q- To do this, define a pointwise inner 
product on tangent spaces TgQ by 

7Q{Vg,Vq):= min Aig{^,0} ■ (2-19) 

We refer to [You 10] for a rigorous treatment of the infinite dimensional case of diffeomor- 
phism groups. We define the vertical subspace of Q at q as 

= ^ dl^Qiq) = 0} (2.20) 

and the horizontal subspace as the orthogonal complement = [q^]^ ■ Denote the orthog- 
onal projection onto 0^ by ,^ 1— )■ {^). This projection depends smoothly on q E Q. The 
vertical projection is similarly written as ^ 1— Vg {^). Let z/i, . . . z/^. be an orthonormal basis 
of . For Vg in TgQ and ^ any generator of Vg, i.e. C,Q{q) = Vg, we can write 



7qK, Vg) = min {7^ (H, (0 + AV,, (0 + AV,) } = min <^ 7g (H, (0 , (0) + V(AT 

= 7G(H,(e),H,(0). (2.21) 

The pair ((5,7q) with 7q defined pointwise by (2.19) is therefore a Riemannian manifold. 
The metric 7q is usually called a normal metric or a projected metric. It coincides with the 
normal metric considered in [GBRll]. 

Tangent action of the vertical space. From (2.16), it follows that Trq^^TQivg)) = Og, for 
^-ll i ^ 0^. Let V^^TQ := ker(T„_^rQ) be the vertical space at Vg. Then irqivg) G Vu^(TQ) can 
be identified with an element of TgQ in the standard way: r : V^^iTQ) 3 ^TQ{vq) ^-> G TgQ, 
where Wg is the unique vector satisfying ^TQ{vg) = ■^\^_q {vq + £uiq)- Let 7q be a Riemannian 
metric on Q and define the Connector of 7q to be the intrinsic map K : T{TQ) — )■ given 
in coordinates as 

i^'ioc(a;, w, u, v) := [x, v + T{x){w, u)), (2.22) 

where F are the Christoffel symbols of the metric. More information on the connector can 
be found, for example, in [MicOS, §13.8]. Its key property is that for any two vector fields 

x,y eX(Q), 

VyX = KoTXoY. 

Moreover, its restriction to vertical spaces Vy^ C T{TQ) corresponds to r, as Kioc{x, e, 0, v) = 
(x, v) = Tioc{x, e, 0, v). For ^ G 0^ one therefore obtains 

r i^TQivq)) = K i^TQivq)) = K {TiQ{q){vg)) = V.^^Q- (2.23) 
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Riemannian submersion. For {G,'yG) and {Q-,'^q) as above, fix a G Q, and consider tlie 
principal bundle projection 

Ii:G-^Q, g^ga. (2.24) 
For all g & G, we decompose the tangent space TgG into the vertical space T^G := kerTgll 

and the horizontal space T^G := (TJ^G)"*" determined by the Riemannian metric 7^. These 
spaces are translations of appropriate subspaces of q, 

T^G = TLgiQ:)=TR,iQ^), and T^^^G = TR^igf) , (2.25) 

where q = Il{g). The second equality in (2.25) is due to the right-invariance of 7^. This 
justifies the terminology vertical and horizontal subspaces for the isotropy subalgebra at 
q and its orthogonal complement . We also record that 

TgUi^cig)) = ^Q^g)) , for all ^7 e G, e e . (2.26) 

It is well known that 11 : G — t- Q is a Riemannian submersion, i.e., 11 is a surjective submersion 
and for any Vg,Wg G T^G, we have 

laKWg) = 7Q iTgU{Vg),TgU{Wg)) . (2.27) 

This property will be useful when we compute covariant derivatives for normal metrics in 
the next section. 



3 Covariant derivatives 

The main goal of this section is to obtain expressions for the covariant derivative of a curve 
q{t) e Q, where Q is equipped with a normal metric. The strategy is the following. First 
we compute covariant derivatives of curves in Lie groups with right-invariant metrics. Then 
we exploit the fact that the projection mapping 11 introduced in Section 2.3 is a Riemannian 
submersion. 



3.1 Preliminary remarks 

Let ((5,7q) be a Riemannian manifold, and let q{t) be a curve in Q. On an open set U G Q 
introduce a coordinate map (j) : U 3 q ^ x E M"' and let V{t) = f*(t) -^\^^^-^ G Tq(t)Q be a 
vector field along q{t). The covariant derivative, with respect to the Levi-Civita connection, 
of V along q{t) is 



Dt y ' ^3 ' Q^k 



(3.1) 



where the Ff^- are the Christoffel symbols of 7q. The geodesic equation -^g = becomes, in 
coordinates, 

+ F*;.x*x^' = 0. (3.2) 
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Let g{t) G G be a curve in a Lie group (Gj'Jg) with right-invariant Riemannian metric 'Jg- 
The symmetry reduced form of the geodesic equation is the Euler-Poincare equation [MR03] 

e + adje = 0, with g = ^Gig), (3.3) 

where ad^ is the metric adjoint, i.e., it has the expression 

ad^K := (ad*(/t^))» (3.4) 

for any z/, k G fl. 

We also recall a formula for the covariant derivative of horizontal vector fields for Rie- 
mannian submersions, which we will use subsequently. Let LI : [Q, 7q) — )■ {Q, jq) be a 
Riemannian submersion and denote the covariant derivatives with respect to the Levi-Civita 
connections on Q and Q by V and V, respectively. Let X,Y ^ ^{Q) be the horizontal lifts 
of X, y G X((5), respectively. Then (see e.g., [Lee97]), 

V^y = V^+^[X,F]^, (3.5) 

where the superscript V denotes the vertical part. The horizontal lifting property of geodesies 
follows. Namely, if q{t) G Q is the horizontal lift of a geodesic q{t) G Q, that is, V^g = 0, 
then q{t) is a geodesic, since V^g = V^g = 0. Note that applying TLI to both sides of (3.5) 
gives 

Tn(v^F)=VxF. (3.6) 
3.2 Covariant derivatives for normal metrics 

The following proposition is a compilation of well-known expressions that will be used ex- 
tensively in the rest of the paper. The proofs of (3.7) - (3.9) can be found, for example, in 
[KM97]. We note that the expression (3.12) below for the horizontal generator was also used 
in [VTIO] and [GBR 11]. 

Proposition 3.1. Let {G,'yG) be a Lie group with right-invariant metric, acting transitively 
from the left on a manifold {Qi'^q) with normal metric 7q. 

(i) Let g{t) he a curve in G, and define ^(t) G by g = ^g{9)- Then, 

(e + adje)^((7). (3.7) 

(ii) More generally, let V{t) G Tg(^t)G be a vector field along a curve g{t) G G. Define 
curves ^(t), v{t) G by g = ^cio) md V = VGid), respectively. Then, 

= + ^ adt + i adU - i^])^{9). (3.8) 

Furthermore, let M{t) G T*^^-^G be a covector field along g{t) and define ^{t) G Q* by 
12 = {TRg)*M. Then, 

= (TV0*(A- ^(ad^(/.«))^ + l^d;,e + ^ad^/i). (3.9) 
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(iii) Let q{t) be a curve in Q and let ^(t) G be a curve satisfying q = CqIq)- Then 

= + adi^(^) H, (O) (g) + V,( V, (0 )q, 



Dt 

In particular, if C,{t) G g^^-j is the unique horizontal generator of q(t), then 

D 
W 



(3.10) 



(3.11) 



(iv) Let q{t) be a curve in Q. The unique horizontal generator of q{t) is given by the Lie 
algebra element J{q) defined by 



3{q):={3{q')y EQf, (3.12) 

where J is the cotangent lift momentum map defined in Section 2.3. In particular, 

D 
Dt' 



(9,J(g) + adt(.^J(g)) (g). 



(3.13) 



Proof. We refer to [KM97] for the proof of (i) and (ii). In order to show (iii) recall the 
projection mapping U : G ^ Q, g ^ ga, for a fixed a E Q. Let q{t) be a curve in Q and 
define the curve ^(t) G to be its horizontal generator, that is, q = Cq{q) and ^ G 0^. Choose 
go G n~^(g(0)) and define g{t) E G by g{fS) = go and g = ^cig)- Then g{t) is the horizontal 
lift of q{t) through go. We apply Til to (3.5) and use (3.7) and (2.26) to find 

D 



q = V^q = TgU {Vgg]= TgU (^ + ad Odg) = + Oq(?)- 



Dt^ 

This shows (3.11). 

Consequently, we have 
D 

Dt' 

so it remains to show that idtVq^^))^ (g) = — Vq( Vg (^) )q. 

We first show that for ri{t) G 0^^^ and with g = g(0), we have 



(3.14) 



q = (dtH.iO + adlj^(^) H, (O)^ (g) = (dt^ + adlj^(^) H, (O)^ (g) - (5* V,(0)q (g), 



d 
dt 



D 
Dt 



t=o 



m)Q{q)- 



(3.15) 



Writing r]{t) = ri{0) + z^(t), we have ?7(0)Q(g) = 0, z/(0) = and we compute 



D 
Dt 



t=0 



D 
Dt 
D 
Dt 

K 
K 



+ ^(t))Q (q) 



i=0 

t=0 
d 

dt 
d_ 
de 



D 
Dt 



ivm)o (q) + (K^))o {q) 



t=0 



(z/(t))o (g) = K 



d 
dt 



Ht))o {q) 



t=0 



t=0 



d 

de 



exp{eu{t)) 



(q) 



e=Q 



K 



d 
de 



e=0 



em)Q (q) = im)Q (q) 



e=0 



d 

dt 



cxp{eiy{t)) 



(q) 



t=0 



vit)] (q): 

t=0 / Q 
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where in the seventh equahty we used that i/(0) = 0. This demonstrates equation (3.15). 
We now show that 



D 

Dt 



t=0 



D 

De 



ivit))a{q{t + e)) = -V^VQ, for r/(t) G gj,) (3.16) 



e=0 



We have 



e=0 



d 

de 
d 

d£ e=0 L 



d 



exp{er]{t)) 



dt 



exp(£r?(t)) 



I^TTQ 



d d 
dt de 



d 
ds 



s=0 



e=0 



ds 



d_ 
de 



exp(er;(t+s)) 



m) 



£=0 



v{s + t)Q{q{t)) 



s=0 



-kttq [ Vero^,,) 



D 
Ds 



vis + t)Q{q{t)) 



s=0 



- Ver^ 



m 



D 

Ds 



V{s + t)Q{q{t)), 



s=0 



where kttq '■ TTQ — )■ TTQ denotes the canonical involution, which in local coordinates is 
expressed as 

K-TTQ{q, g, Sq, 5q) = (g, 5q, q, 5q) . 

Since K {v{t)TQ{qm = K {Tr]{t)Q{qm = WeLo (^(^ + ^)) ^ ° ^er,, = tdrg, 

we obtain the desired formula in (3.16) above. 

Combining the results in (3.15) and (3.16), we get (<9tV^(^))g (g) = — Vg( Vg (^) )^ which 
proves claim (iii) of the proposition. 

The final step is to prove (iv). For a fixed g G Q, arbitrary Wg, Vg G TgQ and ^ G with 
^Q(g) — '^q obtain the following chain of equalities; note that (J(w^))'' G Qg is horizontal 
because of (2.18). 



7Q(((JK))«)Q(g),^,) =7g((J(K))^0 = (JK),0 

= W'eQ(g)>^.Q,^Q 



1' ' '/ g*xg 

lQ{Wg,Vq). 



(3.17) 

(3.18) 



Since Vg was arbitrary we conclude Wg = {{3{Wg)Y)Q{q). This, together with (3.11), shows 
(iv). ■ 



4 Cubics for normal metrics 

In this section we derive the equations of Riemannian cubics for normal metrics. For ease 
of exposition we first consider a more general context and then particularize to the case of 
Riemannian cubics. The examples of Lie groups and of symmetric spaces are worked out in 
detail. 
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4.1 Preparations 



Consider a manifold Q with a linear connection on its tangent bundle TQ. Denote the 
covariant derivative with respect to this linear connection by For Vg, Wg G TgQ write 
(Vg)^ G T^qTQ for the horizontal lift of Vg to Wg, i.e., in local coordinates, 



where F is the Christoffel map of the linear connection. The vertical lift of Vg to Wg is written 
^^(t, s) G Tgi^t){TQ) splits into horizontal and vertical parts 

5q = {6q)fg^g^ + ( 



{'^q}w '■— ^"^q)' "^^^ ^ Variation (t, s) i— )■ q(t, s) of a curve q{t) = q{t,0) the curve 



(4.1) 



For a function ^ : TQ — t- g and an arbitrary Vg G TQ define the g- valued linear form |^|^ by 

(4.2) 



<5g 







de 




A 




de 


L=o 



^iv{e)), for any G TgQ, 



vie 



[w. 



^ . On the other hand we write 14 1 



where v{e) is any curve in TQ with ^^i^^g-v-y - v-g;^,,• — — — 
for the fiber derivative of ^ at Vg. Note that if q{t, s) is a variation of a curve g(t) = q{t, 0), 
then using the splitting (4.1) we get 



s=0 



(4.3) 



We define the operator {^1^)* '■ d* ^ T*Q by 



6q 



fl,Wg 



fx, 



5g 



,UJg 



, for any G fl*, G T^Q, (4.4) 



and similarly the operator (^|^ )* • S* ~^ T*Q by (keeping track of q and q) 



6q 



fl,Wg 



5g 



for any /x G 0*, Wg G T^Q. (4.5) 



4.2 A generalized variational problem 

Given a Lagrangian £ : 20 — t- M and a smooth map ^ : TQ — 0, consider the action functional 
on the space of curves q{t) : [0, 1] ^ Q given by 



J{q] 



(4.6) 
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and Hamilton's principle SJ^ = with respect to variations satisfying Sq{0) = Sq{l) = and 
Sq{0) = Sq{l) = 0. As we shall see below, the cubic spline Lagrangian for normal metrics fits 
in this framework. 

Taking variations of J' we obtain, using (4.3), 



5J 



5e , A /6£ , , 

6£ _ d 6£ /5i 
5^ ~ 'dtJl' \Tq 
5i\* f5£ _ d 6i 

Jq) \M~JtTi 

6qJ Dt \6q 



Si d Si 



dtS^ 
S^ D 



S^ ) dt 



Sq Ds 



q 



s=Q 

Sq 

Si _ d Si 
si " diJi 



dt 

Si _ d Si 

si ~ ~dtTi 



D 
Wt 



Sq ) dt 



,Sq) dt 



The Euler-Lagrange equation for Hamilton's principle SJ^ = is therefore 

E d S 



D Si 

o 

Dt \Sq 



0. 



4.3 Cubics for normal metrics: Euler-Lagrange equations 



(4.7) 



Let (G, '-)g) be a Lie group with right-invariant metric 7^, acting transitively from the left on 
a manifold (Q, 7q) with normal metric 7q. Let q{t) be a curve in Q, originating at go = ^(O). 
Recall from (3.12) that its horizontal generator curve is given by J(g) = (J(g''))'' G Q^. 

Lemma 4.1. For any curve q{t) G Q, the curve dt{^{q)) + adj^^^ J(g) is horizontal, that is, 

Proof. For any Qq G n~^(go) define the horizontal curve g{t) G G by g{0) = go and g = 
(J(g))^ (g). Then, by formula (3.5), ^g is horizontal. Moreover, by Proposition 3.1, 



D 

Dt' 



-g = [d,3iq) + adt^^^ J(g))^ (^?) = TR, (9,J(g) + adt^.^ J(g) 
The statement of the lemma now follows from (2.25). 



This lemma enables us to rewrite the Lagrangian (2.8) of Riemannian cubics, evaluated 
along the curve q{t), as follows. 



L{qA,<i) 



1 


D 


' 1 


2 


Di^ 


."2 



d,3{q) + e.dl^^^3{q)j^{q) 



9,J(g) + ad^(^)J(g) 



Hamilton's principle (2.5) for Riemannian cubics is SJ^ = with cost functional J7 of the 
form (4.6), 



J[q] 





D 




f- 


Di"^ 


dt= [ 


/o 2 


q Jo 



i[m.dtm) dt, 
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where 



6 + adl 6 



Remarkably, the function i coincides with the reduced Lagrangian in the Euler-Poincare 
reduction of Riemannian cubics on Lie groups, see (4.17) below. The variational derivatives 
of i read 



= (adgi vf - ad; and ^ 

Hence, the Euler-Lagrange equation (4.7) becomes 
D 



where rj := C,2 + ad^^ 



Dt 



63 

Tq 



{-dtVi + (adj r]f - ad; J^) = 0, where := J + ad^ J. (4.8) 



We emphasize that the covariant derivative is understood with respect to a chosen linear 
connection on the tangent bundle TQ. A possible choice is the Levi-Civita connection with 
respect to the normal metric 7q on Q. 



4.4 Splines on Lie groups 

Riemannian cubics on Lie groups (G, 7g) with right-invariant metrics 7g were treated in 
[GBHM+10] by second-order Euler-Poincare reduction. Here we revisit the problem from 
the point of view of normal metrics. 



Equations of motion. We first observe that 7^ is a particular case of a normal metric. 
Namely, let (G, 7^) act on G by left multiplication. Then the normal metric induced on G 
is again 7^. The generator of a curve g{t) is the right-invariant velocity vector, 3{g,g) = 
TRg-ig. Choosing as linear connection on TG the Levi-Civita connection of 7g one arrives, 
by Proposition 3.1 (ii), at 



53 



i9,g) 



: Sg ^ -^adl^ -^adlri -hc,7]], where ^ ■= TRg-ig, v := TRg-iSg, 



5J 



:/i^^(ri?,-i)*(^^ad;«e'-^(ad5/i»)'-^ad*/i^, where ^:=TRg^ig, 



63 

Tg 



: Vg TRg-iVg for any Vg G TgG, and 



{9,g) 



53 



-.ft^iTRg-iYfi. 



{9,9), 



Using now the expression for given in (3.9) we get 



63 

Tg 



D 
Dt 



53 



^i={TRg-iy{-dt-ad*j)^i 



(4.9) 
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for any curves g(t) G G and fi{t) G Q*. The Euler-Lagrange equations (4.8) are therefore 

{TRg-i)*{dt + ad*j)[dtV^ -{adj7]f + a.d*T] = 0, where r/ := j + ad^ J. (4.10) 
This is equivalent to 

(9t + adj) [Str^^- (adjr7)^ + ad;j^] =0, where r]:=j + adtj, (4.11) 
and J := 3{g,g) = gg-\ 

Bi-invariance and the NHP equation. If the metric 7^ is bi- invariant, then ad^ = — ad, 
and the above expressions simphfy. Namely, 77 = J and (4.11) is equivalent to the NHP 
equation 

J + [J,J] = 0. (4.12) 

This equation was originally derived in [NHP89] for the special case of G = S0{3) with 
bi-invariant metric. It can be integrated once to yield 

J + [J,J] = z/ (4.13) 

for a constant u E g. Solutions J were called Lie quadratics in [Noa03b, Noa04]. For 
G = 5*0(3), long-term behavior and internal symmetries of Lie quadratics were studied 
there, both in the null [u = 0) and the non-null (z/ 7^ 0) case. Generalizations to cubics in 
tension can be found in [NP05]. 



Euler— Poincare reduction. We include a discussion of second-order Euler-Poincare re- 
duction and, in particular, the Riemannian cubics in this context. More details and higher- 
order generalizations can be found in [GBHM+10]. The first-order case is discussed, for 
example, in [MR03]. Start with a right-invariant Lagrangian L : T(2)G R with reduced 
Lagrangian £ : 20 — )■ M. Consider Hamilton's principle 

6 [\{g,g,g)dt = (4.14) 
Jo 

with respect to variations of curves g{t) : [0, 1] — ?■ G respecting boundary conditions 5g{0) = 
Sg{l) = and ^^'(0) = Sg{l) = 0. The right-invariance of L leads to the equivalent reduced 
formulation ^ 

5 / i{^,Odt = 0, (4.15) 
Jo 

with respect to constrained variations of curves ^(t) : [0,1] — )■ Q. Namely, one considers 
variations of the form 5^ = rj] with r]{t) : [0, 1] — ?■ arbitrary up to boundary conditions 

7/(0) = r]{l) = and r]{0) = 77(1) = 0. Solutions g{t) of (4.14) and solutions ^(t) of (4.15) 
are equivalent through the reconstruction relation ^ = TRg-ig = 3{g,g). 
Taking constrained variations of (4.15) leads to the second-order Euler-Poincare equation 

— + adf T- - ^ — r =0. (4.16) 
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For a curve g{t) G G with right-invariant velocity vector ^{t) = TRg^ig we rewrite the 
Lagrangian of cubics, 



Lig,g,g) 



1 


D 


2 


1 


2 




9 


2 



e + adje 



where in the second equahty we used (3.7) and the third equahty follows from right-invariance 
of 7g. This demonstrates that the spline Lagrangian L only depends on the right-invariant 
velocity ^ and its time-derivative ^ and is therefore right-invariant. The reduced Lagrangian 
can be read off as 

1 • +2 



(4.17) 



The dynamics are governed by the second-order Euler-Poincare equation (4.16), which be- 
comes, for i as above, 



{dt + ad^) [dtT]' - (adg r])' + ad; C] = 0, where := ^ + adj 



This coincides with (4.11), since ^ = TRg-ig = 3{g, 



4.5 Splines on symmetric spaces 

We particularize the equation for Riemannian cubics (4.8) to symmetric spaces. Due to the 
appearance of the horizontal generator, this equation lends itself to the analysis of horizontal 
hfting properties to be addressed in Section 5. We also comment on how it is related to the 
equation derived in [CS95]. 



The horizontal generator. Recall that for any curve g{t) in a Lie group G with g = TRg^ 
for ^(t) G 5 and any curve v{t) G fl. 



d 
dt 



Adg-i z/ = Adg-i (z>+ [z/,e]). 



(4.18) 



Let G be a Lie group with a bi-invariant metric 7 that acts transitively on a manifold Q 
equipped with the normal metric 7q, so that the action is by isometrics. Denote by Ga the 
isotropy subgroup of a fixed element a G Q, so that Q is diffeomorphic to G/Ga-, the quotient 
being taken with respect to the right-action of Ga on G. Recall the Riemannian submersion 
n : G — )■ Q given hy g ga. Note that for any g & G with ga = q, 



Adg o Ha o Ad^ 



-1 



and Adg o Va o Ad 



-1 



(4.19) 



Lemma 4.2. Let g{t) E G he a horizontal lift of a curve q[t) G Q, and let J := J(g) he the 
horizontal generator of q{t). Then, 



V, (J) = 0, V, ( J ) = 0, Vg (J + 



J, J 



0, 



V, J + 2 



J, J 



+ 



J, J 



(4.20) 
(4.21) 
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Proof. The first equation follows from horizontality of J. It is equivalent to (Ad^-i j) = 0. 
The second equation follows from taking a time derivative of this last relation using (4.18). 
Note that g = TRgJ. For instance, 

dt {Va (Ad,-i J)) = V, (Ad,-i j) = 0. 

Therefore, Vg ^jj = 0. The third and fourth equations follow from taking two more time 
derivatives. ■ 



Symmetric spaces. Assume in addition that a is an involutive Lie algebra automorphism 
of Q such that and are, respectively, the +1 and —1 eigenspaces. Then {G,Q,a) is a 
symmetric space structure. The following inclusions hold for all q & Q, 



Lflg > 05 J C 0q , [0g , 0g J C 0, , [0g , 0g J C 0g . 



(4.22) 



The first identity follows because is a Lie subalgebra of q. The second one is a consequence 
of the Ad-invariance of the inner product on g. The third one is characteristic of symmetric 
spaces. It follows from the eigenspace structure of a described above. We now compute the 
the Euler-Lagrange equation of cubics, (4.8). We will find that it is equivalent to 



H, ( J + 2 



J, J 



0. 



(4.23) 



We start with a lemma. 

Lemma 4.3. Relative to the Levi-Civita connection on TQ, 



(<?,g) 



.5q ) = [J{6q),j{q)] , and 



3{5q). 



Proof. The second statement is due to the linearity of J on fibers of TQ. In order to prove 
the first equation, let q{e) be a curve with g(0) = q and de=oq{s) = 6q, and let g{e) with 
5^(0) = ghe horizontal above q{6). We construct the parallel transport of q along q{e). Define 

u{e) := J(g) + e [3{6q), J(g)] , X{e) := TRg^,)Uj{e). (4.24) 

To first order in e, X{e) is a horizontal vector field along g{£) and uj{e) lies in 0^£)- Now define 
a vector field along q{e) by v{e) := Tg(^^)Il{X[e)). Note that to first order in e, 3{v{e)) = uj{e). 
Denoting by V the covariant derivative on G with respect to the Levi-Civita connection of 
7, we use (3.8) to get 

1 



l9iSq),3iq 



(9) 



V 



G 



where we used (3.5) in the second step. Recall that [j(5g), J(g)] is in , so that applying 
Til to the above shows VsqV = 0. To first order in e, v{e) is therefore the parallel transport 
of q along q{e), and 



,5q 



d 

de 



3{v{e)) 



e=0 



d 

de 



uj{e)= [J((5g),J(g 



e=0 
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It follows from this and from Lemma 4.3 that 



63 

Tq 



53 

Tq 



^-"^ ((/^%(^)) ' (4-25) 



for any fi G 0*. Taking the sharp of (4.8) therefore gives 



J, J 



D 
Dt 



J J?) 



0. 



(4.26) 



It follows from (4.20) and relations (4.22) that J + 
and therefore 



J, J 



is the horizontal generator of (j) (g) 



D 
Dt 



J J?) 



J + 



D 

' Dt 

J, J 



J + 

1 
2 



J, J 



(?) 



J, 



J, J 



(?) 



J + 



J, J 



Q 



(?), 



where in the last equality we used (4.22), namely 

([J, [J,3]])^{q) = -([J, V, (J)])^(g) = -([J,5])^(g). 
Therefore, (4.26) becomes (4.23), 



H, ( J + 2 



J, J 



0. 



(4.27) 



(4.28) 



We will exploit the similarity of this equation with the NHP equation (4.12) when we analyze 
the horizontal lifts of cubics in the next section. 

Remark 4.4. The equations for Riemannian cubics in symmetric spaces were first given in 
[Cyy-j]. We briefly remark on how those equations are related to (4.28). We first note that 



(4.21) together with (4.22) show that the vertical part of J + 2 



J, J 



vanishes. Hence, 



3+2 



3,3 



Let g{t) be a horizontal lift of a Riemannian cubic g(t), and define V(t) G by 



V = Adg-i J. 



One checks easily that (4.29) implies 



V + 



V, 



v,v 



(4.29) 



(4.30) 



(4.31) 



This coincides with equation (46) of [CS95]. 
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Example: G = S0{3) and Q = S'^. Let G = S0{3) and Q = CR^ and let ^0(3) act 
on through its action on vectors in M^. We start with our notational conventions. 

Remark 4.5. [Conventions for 5*0(3) and S'^] 

Throughout this paper we use vector notation for the Lie algebra so (3) of the Lie group of 
rotations S0{3), as well as for its dual 5o(3)*. One identifies so(3) with via the familiar 
isomorphism 



^so(3), ri= 6^fi:=ri= a -c , (4.32) 





called the hat map. This is a Lie algebra isomorphism when the vector cross product x is 
used as the Lie bracket operation on M^. The identification of so(3) with induces an 
isomorphism of the dual spaces so (3)* = (M^)* = M'^. We represent tangent and cotangent 
spaces of 5*^ as 

T^S^ = {(x,v) G 52 X I X. V = 0} , T;52 = {(x,p)G52xM3|x.p = 0} (4.33) 

with duality pairing ((x, p) , (x, v))^«^2xT52 = P • v. Whenever admissible, we will drop the 
explicit mention of the base point x in what follows. 

The infinitesimal generator is given by (17)52 (x) = f2 x x. We consider the bi-invariant 
extension jsois) to 5*0(3) of the identity moment of inertia inner product (f2, f2)^^j.3^ = 
O • O = ||r2|p on so(3). The corresponding normal metric on 5*^ is the round metric. The 
vertical and horizontal spaces are, see (2.20), 

so(3)^ = {ri Gso(3)|ri = Ax for A G M} and so(3)^ = {fl G so(3) |ri ■ x = O} . (4.34) 

The map J of (3.12) becomes J(x, v) = x x v. Equation (4.28) is 

X X ( j + 2J X J) = 0, with J = X X X. (4.35) 

Equation (4.35) appears in [Kra05], where it is derived from the general Euler-Lagrange 
equation for cubics (2.10). The similarity of (4.35) with the NHP equation (4.12) on 5*0(3), 
J + J X J = 0, was remarked upon there. We will take advantage of this similarity in Section 
5 for the investigation of horizontal lifts of cubics. 



4.6 Curvature from cubics 

One can infer the Christoffel symbols for a given metric if one knows the geodesic equation. 
It is interesting to note that in a similar way one obtains expressions for the (sectional) 
curvature from the equation of Riemannian cubics. We illustrate this in the case of Lie 
groups and symmetric spaces. 
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Lie groups. Consider a Lie group G witli metric 7^ tliat we assume, for simplicity, to be 
bi-invariant. Analogous arguments apply in the non-bi-invariant case. Let g{t) be a curve in 
G and J = TRg-ig its right-invariant velocity. Recall that ad^ = — ad and from (3.7), (3.8) 
that 

D,g = {dtJ)^ (g), Dig = (J - ^[J, 3]) ^{g) =: (z/)a(^7), D^g = [v - ^[J, v])^{g), 

where we defined v := J — |[J, J]. The general Euler-Lagrange equation for cubics (2.10) 
becomes 

u-\[3.u]^^{g) + R{3G{g)JG{9)) (Jg(^7)) =0. (4.36) 
On the other hand, if g{t) is a spline, then the NHP equation (4.12) is satisfied, and therefore 



J, J 



. Plugging this into (4.36), yields 



J, J 



G 



(^?) + i?(jG(^),JG(^7)) (Jg(^7)) =0. 



We conclude that for any g E G and C,,i] G Q, 



R{vg{9),^g{9))^g{9) = -^([e, [^,v]])ai9)- 



(4.37) 



Symmetric spaces. For symmetric spaces one derives in a similar fashion that for a cubic 
q{t) with horizontal generator J = J{q), 



DU = J + 



J, J 



It follows from (2.10) and (4.28) that 



R{Dtq,q)q 



J, J 



J, J 



We conclude that for any q E Q and 77, ^ ^ Qq 1 



R{m{l).Ul))iQ{l) = -mi.r,]])Q{q). 



5 Horizontal lifts 



The horizontal lifting property of geodesies in the normal metrics context has been an im- 
portant feature of the large deformation matching framework in computational anatomy. 
References [HRTY04, MTY06, YAM09], amongst others, explore this aspect in detail. This 
motivates us to ask which Riemannian cubics on the manifold Q lift to horizontal cubics on 
the Lie group G. In the case of symmetric spaces we give a complete characterization of the 
cubics that can be lifted horizontally. We then consider the more general setting of Rieman- 
nian submersions and formulate necessary and sufficient conditions under which horizontal 
lifts are possible. 
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5.1 Symmetric spaces 



Here we completely characterize the Riemannian cubics that can be lifted horizontally to the 
group of isometrics. Let {G, Q, a) be a symmetric space structure, as defined in Section 4.5. 
In particular we recall the important relations (4.22), 



[0g , 0g J C 0, , [0g , 0g J C 0^ , [0g , 0g J C 



(5.1) 



Theorem 5.1. A curve q{t) & Q is a Riemannian cubic and can he lifted horizontally to a 
Riemannian cubic g(t) E G if and only if it satisfies q{t) = (C(^))q(Q'(^)) for a curve ^(t) G 
of the form 



ut' 



+ vt + w, 



(5.2) 



where u,v,w span an Abelian subalgebra that lies in 0^o)- 



Proof. Suppose q{t) is a cubic and can be lifted horizontally to a cubic g{t). The horizon- 
tal generator curve J(g(t)) simultaneously satisfies equations (4.12) and (4.29). Therefore, 
[J, J] = 0. In particular 

2 



7 



J, J 



7 J 



J, J 



J, J 



where we used (4.20) and (5.1). We conclude that [J, J] 
equation (4.13) this reveals that J is constant. Therefore, 



0, 



0. Together with the NHP 



J(g(t)) 



uV 



+ vt-\- w 



with constants u,v,w G 0. These are mutually commuting because of [J, J] = and its 
time derivative [J, J] = 0. Moreover, J and J are horizontal. Due to [J, J] = 0, J is also 
horizontal. Therefore, u,v,w E 0^o)- Setting ^{t) = 3{q{t)) therefore completes the first part 
of the proof. For the reverse, let q{t) be a curve that satisfies q{t) = {C,{t))Q{q{t)) for C,{t) G 
of the form (5.2) with mutually commuting u,v,w G 0^o)- element a E Q and recall 

the Riemannian submersion 11 : G — )■ Q given hj g ga. Assume without loss of generality 
that g(0) = a. We first show that ^{t) is horizontal at all times. Then ^(t) and ^(t) are 



also horizontal, since the commutators 



and 



vanish. Define the curve g(t) G G 



by g{0) = e and g = TRg^. This curve lies in the Abelian subgroup Exp(span(-u, f , w)) of 
G with Lie algebra span(M,f, w) C 0^o)- Therefore, TLg-ig = Ad^-i ^ G span(-u, t>, w) which 
means that ^ G Ad^ span(-u, f , w) C Adg0^Q^ = 0^j)- Hence ^(t) is the horizontal generator 
of q{t), and therefore g{t) is a horizontal lift of q{t). Moreover, g{t) is a cubic in G since C,{t) 
satisfies the NHP equation (4.12). It also satisfies (4.28). Therefore, q{t) is a cubic in Q. 
This concludes the proof. ■ 

A Cartan subalgebra ( CSA ) based at q E Q is a maximal Abelian Lie subalgebra contained 
in 0^. The rank of the symmetric space is the dimension of its CSAs. The greater the rank 
of a symmetric space, the larger the set of vectors u, v, w consistent with the requirements of 
Theorem 5.1. 
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Corollary 5.2. In rank-one symmetric spaces the only Riemannian cubics that can be lifted 
horizontally to Riemannian cubics on the group of isometrics are geodesies composed with a 
cubic polynomial in time. 

Proof. Let Q be a rank-one symmetric space. Then any CSA is one-dimensional. A curve 
q{t) is therefore a cubic that can be lifted horizontally if and only if q{t) = {^{t))Q{q{t)) with 
^(t) €0 of the form 

m= (^"^ + bt + c^d, (5.3) 
where d G g^Q-, and a,b,c G M. Therefore, 

q{t) = e(^+'^+^*)'^g(0), 

which corresponds to the geodesic y{t) = e*'^g(0), composed with a cubic polynomial in 
time. ■ 

Remark 5.3. Consider G = 5*0(3) and the rank-one symmetric space Q = S'^. The special 
cubics appearing in Corollary 5.2 are considered in more detail in [PR97]. In particular, the 
case b = c = corresponds to the so-called natural splines in computer aided design (CAD) 
applications. 



5.2 Riemannian submersions 

In this section we generalize the question of horizontal lifts of cubics to the Riemannian 
submersion setting. We also show how the result implies Theorem 5.1 of the previous section. 
The result is the following: 

Theorem 5.4. Let U : Q ^ Q be a Riemannian submersion, and let q{t) E Q be a Rieman- 
nian cubic. Moreover, let q{t) E Q be a horizontal lift of q{t). The curve q{t) is a Riemannian 
cubic if and only if 

D 



'^Dt' 



= 0, for all t, (5.4) 
where the superscript V denotes the vertical part. 

Remark 5.5. For any two horizontal vectors v,w E T?Q, the expression is defined 

as [-D,?!;]^ := [V^,VF]^(g), for horizontal extensions V ,W of v.,w. 

Proof. We denote v{t) := q{t) and v{t) := q{t). The metrics on Q and Q are denoted 7 and 
7. The covariant derivatives with respect to the Levi-Civita connections are written V and 
V respectively. Recall that by definition = V„f and -^v = VyV. Using the formula 

= + \ [X, Y]^ for the covariant derivative induced by a Riemannian submersion 

of horizontally lifted vector fields X and Y, one obtains 

= v7v, i^.yv = {y^v + hy, ^vY, (5.5) 

= + (^]V + 1 v,([5, V^vf). (5.6) 
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Suppose q{t) and q{t) are as in the statement of the theorem and let q(t) be a Riemannian 
cubic sphne, i.e., the respective sphne equations are satisfied, 



{V^fv + R{V^v,v)v = 0, 



(5.7) 



For the following manipulations we record that 



d 
dt 



where the second step follows since VyV is horizontal. We use this equality as well as equations 
(5.5) - (5.6) to obtain 



1. 



7((Vc)^C,Vc5) =^(^{V,yv,V,vj +-7 (V^i[v,V,vY),V,v 
= 7 ((V^t^, V::^;) - ^7 ([^^, V:^]'', (V,)^^ 



7((V,)3t;,V,t; 



[V, VyV] 



Hence, 



7 (Vc)^^, Vci) = 7 {C^vfv, V^v) - - [v, 



{5.i 



On the other hand O'Neill's formula for sectional curvatures of Riemannian submersions 
[0'N66] (Equation 3. in Corollary 1), implies that 



[v, VvV] 



(5.9) 



7 \^R{VvV, v)v, VvVj = 7 {R{VyV, v)v, Vvv) - - 
Adding (5.8) and (5.9) and using the spline equations (5.7) we conclude that 

[v,v7vV = [i,^,v]'' = 0, (5.10) 

which is (5.4). To show the reverse direction, we note that when (5.4) holds, then (5.5) and 
(5.6) take the simplified form 



2~lV 



(5.11) 



Hence, (Vc)^w is horizontal. Moreover, (Vi})^^ splits naturally into horizontal and vertical 
parts. Therefore, checking that the second equation in (5.7) holds, amounts to verifying that 
for any choice of horizontal vector field h{t) and any choice of vertical vector field w{t) along 



7 (r{VvV, v)v, -1 (Ri^vV, v)v, h) = 
7 (^R{ViV, v)v, + |7 [[v, {ViYv]^, w 



0, 



(5.12) 
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where we denoted h := U^h. To proceed, we introduce the (1, 2)-tensors A and T defined, 
for arbitrary vector fields E, F, by 

AeF = (Veh(F^))'' + {Veh{F'')Y (5.13) 

TeF = (vEviF""))" + (vEv{F'')y . (5.14) 

The superscripts H and V denote the horizontal and vertical parts, respectively. Definitions 
(5.13) and (5.14) coincide with the ones given in [0"N66]. It is shown there (in Equations 
{3} and {4}) that if X, Y, Z, H are horizontal vector fields and W is a. vertical vector field, 
then 

7 (^R{X, Y)Z, if) = 7 {R(Jl*X, UX)'n*Z, U,H) + 27 {AxY, AzH) 

- 7 {AyZ, AxH) - 7 {AzX, AyH) . (5.15) 

and 

7 (i?(X, Y)Z, l^) = - 7 {{yzA)xY, l^) - 7 {AxY, TwZ) 

+ 7 {AyZ, TwX) + 7 {AzX, TwY) . (5.16) 

Note that we differ from [0'N66] in our sign convention for the curvature tensor. It is also 
shown in [0'N66] that for any two horizontal vector fields X and Y one has AxY = ^[X, Y]^ . 
In particular, (5.10) can be written as AyCVvv) = 0. This, together with (5.15), implies that 

7 (j^iVaV, v)v, /i) = 7 {RiyvV, v)v, h) , 

which is equivalent to the first equation in (5.12). In order to show the second equation we 
take a covariant derivative of (5.10) written in the form A:{,{Wvv) = to obtain 

= VMv{^vv)) = (VsA)s(Vc5) + A,{{Vi,fv). 

It follows from this and (5.16) that 

7 (^R{'VvV,v)v,'LL?j = -7 (jk{v,'VvV)v,W^ = 7 (^{'VvA)ii{'VvV),W^ 

= -7 [A,ii\/,yv),w) = -^7 {[v, (V,)^^]^,*) . 
Therefore the second equation of (5.12) is satisfied. This concludes the proof. ■ 

5.2.1 Example: Normal metrics in the bi-invariant case 

We now show how this result relates to Theorem 5.1 of Section 5.1. Let G be a Lie group 
with a bi-invariant metric 7 that acts transitively on a manifold Q equipped with the normal 
metric 7q. Recall that for a fixed element a & Q the map II : G ^ Q , g ^ ga is a, Riemannian 
submersion. 
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Lemma 5.6. Let g & G with Il{g) = q and let ^ and rj be in . Then 



M9UGi9r = iy,i[^,v]))Gi9)- (5.17) 

Proof. In order to compute the left hand side we extend ^cig) to a horizontal vector field 
on G. Namely, H^{h) = /acI _i(5)(^)) where denotes the left-invariant vector field lu{h) = hv 
for any z/ G 0. Similarly H^{h) = /Ad _i(»?)(^)- Now 

T/ / \ V 



D 



MgUcigW = [h^.h^Y (g) = [Ia^^^a^M) = i[^,v]Gig)y = (v,([e,r7])),,(^?). 

■ 

Theorem 5.7. Letq{t) & Q be a Riemannian cubic spline with horizontal generator 3 {q(t)) G 
g^^-j. A horizontal lift g{t) G G of q{t) is a Riemannian cubic if and only if 

Vg ([J, J]) = 0, for allt. (5.18) 

Proof. Since g = (J)^ (g) and DfCj = (3^^{g) with both J and J in , Lemma 5.6 gives 

{y,{[3,3]))jg). (5.19) 

This expression vanishes according to Proposition 5.4, and therefore Vg ^[J, J]^ =0. ■ 

Remark 5.8. If Q is a symmetric space, then relations (5.1) hold. Therefore Vg ^[J, 3]j = 
is equivalent to [J, J] =0. This property was the basis for the proof of Theorem 5.1. 

6 Extended analysis: Reduction by isotropy subgroup 

In the previous section, we analyzed the relationship between Riemannian cubics on Q and 
horizontal curves on G. More precisely, we gave a necessary and sufficient condition guaran- 
teeing the existence of a horizontally lifted cubic on G covering a given cubic on Q. In the 
present section we include in our considerations the non-horizontal curves on G. We show 
that certain non-horizontal geodesies on G project to cubics. We then extend the analysis in 
the following way. We reduce the Riemannian cubic variational problem on G by the isotropy 
subgroup Ga of a point a E Q. The reduced Lagrangian couples horizontal and vertical parts 
of the motion, which accounts for the absence of a general horizontal lifting property. The 
reduced form of the equations reveals the obstruction for a cubic on G to project to a cubic 
on Q. The main technical tool in this section is second-order Lagrange-Poincare reduction. 
The main references are [CMROl] for the first-order theory and [GBHRll] for the recent 
generalization to higher order. 
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6.1 Setting 

Let G be a Lie group with bi-invariant metric 7g, acting transitively from the left on a 
manifold Q with the normal metric 7q. Recall the Riemannian submersion U : G ^ Q given 
hj g ga for a fixed element a E Q. Let Ga C G he the stabilizer of a. Consider the right 
action 

ip : G X Ga ^ G , {g,h) gh. 

This action is free and the projection from G onto the quotient space G/Ga is a submersion. 
Moreover the map G/Ga ~^ Q given by [g] ga is a diffeomorphism. The ingredients 
{G,Q = G/Ga,Ga,^,'ip) therefore constitute a principal fiber bundle with total space G, 
base manifold Q, structure group Ga, projection 11, and action ip. The Lie algebra of the 
structure group Ga is Qa = Qa- Recall the vertical and horizontal projections, Vg '■ Q ^ 
and Hg : — 7- 0^, for any q E Q. Note that for any g E G with ga = q, 

AdgoH^oAd;^ = and Ad^ o o Ad^ ^ = . (6.1) 

We equip G with the principal bundle connection A, 

A:TG^Qa, Vg^ Ag{vg) := V,{TLg-iVg). (6.2) 

Recall from Section 2.3 that 7^, together with the projection 11, induces a splitting of TG 
into horizontal and vertical subbundles TG = T^G © T^G. This splitting coincides with 
the one prescribed by the connection A, that is, T^G = ker^^ and T^G = kerT^II. The 
curvature of A is the g-valued two-form 

Vg) = {Va{TLg-lUg), Va{TLg-lVg)] - Va {[TLg-lUg, TLg-lVg]) 

= -Vai[Ha{TLg-iUg),HaiTLg-iVg)]) ^ • ^ 

for Ug,Vg G TgG. Cousidcr the following action of Ga on G x g^, 

GaX{GxQa)^Gx 0„, {h, g, ^ (gh, Ad^^ 0- (6-4) 

We define the associated adjoint vector bundle over Q, Qa := (G x Qa)/Ga- The equivalence 
class, i.e., orbit, of {g,C,) G G x will be denoted by square brackets, a = [g,C,] G Qa- The 
principal connection A induces a linear connection on Qa with covariant derivative 

^m,m] = [9{t),m + [A{m),m]] ■ (6.5) 

We will sometimes use the shorthand a := where a{t) is a curve in Qa- Moreover, we 
define the map 

i-Qa^Q, [g,7]]^ ■i{[9,v]) ■= -^(^gV, (6.6) 

and write [g,"!]] =: o" i— ?> a := i[[g,ri]), as shorthand. Note that i{[g,ri]) G 0^. We introduce 
the fiber-wise inner product 7 on Qa given by 

7(a,p) := 7(0-, p) 

and its corresponding norm is denoted by ||.||g„. We define the 0^- valued reduced curvature 
2-form B, 

B{ug,Vq) := [g,B{ug,Vg)], 
for Uq,Vg G TgQ, where g E G and Ug,Vg G TgG are such that 11(5') = 1 ^^'^ '^g^i'^g) — ^<?' 

TgU{Vg) = Vg. 



Gay-Balmaz et al. Invariant higher-order variational problems II 



29 



6.2 First-order Lagrange Poincare reduction 

We start by recalling first-order Lagrange-Poincare reduction, which makes use of the bundle 
diffeomorphism 

a'i> : TG/Ga ^ TQ Xq~q,, [g,g] ^ (g, g) x [g,A{g)]. (6.7) 

Here we introduced the quotient TG/Ga of TG by the natural action of Ga, whose elements 
we denote by [g,g] G TG/Ga- We also defined {q,q) := Tgll{g,g) for any representative 
{g,g) of [g,g]- Let L : TG — )■ M be a Ga-invariant Lagrangian. The reduced Lagrangian i 
is defined by L = £ o aj^\ In order to compute the Lagrange-Poincare equations one takes 
constrained variations in the reduced variable space. Take for instance the kinetic energy 
Lagrangian L{g,g) = ^ \\g\fg. Hamilton's principle ^S* = for S* = Ldt leads to the Euler- 
Lagrange equation Dtg = 0. This is the geodesic equation on G. We derive the corresponding 
Lagrange-Poincare equations. The reduced Lagrangian i : TQ XgQa is given by 

= I \\q\\l + I \W\\l = I \\q\\l + I Ml , (6.8) 

where we recall that a := ^(0"). We will need the following result. 

Lemma 6.1. Consider the map i : Qa ^ Q defined in (6.6) and a curve G Qa, covering the 
curve q^ E Q. Then, we have the formula 

where the covariant derivative D'^ was defined in (6.5) and the map F : TQ Xq Qa ~^ Q is 
defined by 

where Vg G TG and rj E Qa are such that TgTC (vg) = Vg and [g, rj] = (Xg. 



De 



In short, 5a = i{6a) + F{6q,a), where 6a := ^ a. Note that i[5a) G 3^ and 

^ e=0 



F{6q,a) G Qg. One computes 

= jQ{q,D^q^) +-f{a,6a)dt = --fQ{Dtq,6q) + 'y{a,6a) dt. 
Jo Jo 



Using the constrained variations 



^(^ = -j^V - W, V] + ISiSq, q) G Qa 



gives the horizontal and vertical Lagrange-Poincare equations 
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where we defined 

7Q {p, S^)" ,Wg^ =1 (P> K^iWq)) = 1 (^P,B{Vg, Wq) 

for all Vq,Wq G TqQ and p G {Qa)q- Using the expression (6.3) one computes that 

(p,i,,^)" = ([JK),p])Q(g) = -V.,pQ, (6.10) 

so that (6.9) becomes 

D 

— g = Vaffo, c^ = 0. (6.11) 

Dt^ ^ Dt ^ ' 

The first of these equations is the geodesic equation on Q up to a forcing term on the right 

hand side. One recognizes again that horizontal geodesies (a = 0) on G project to geodesies 

on Q. 

Remark 6.2. The projections to Q of non-horizontal geodesies in G are called ballistic curves 
in [AKLM03]. 

Remark 6.3. Note that since the Lagrangian L{g,g) = ^\\g\fg is not only Ga-invariant 
but G-invariant, the Lagrange-Poincare equations (6.11) can be further reduced to yield the 
Euler-Poincare equations 



d 

It"--' 



^ = _adU = 0, ^ = TRq-.g. 



We thus obtain that ^ is a constant. This form of the equation is however not adapted for 
our purpose since it does not involve the manifold Q. 

Remark 6.4. If Q is a symmetric space, then the first equation in (6.11) can be used together 
with the second equality of (6.10) to write 

Since both J and [J,o"] are in we conclude that J = — [J,a]. The second equation in 
(6.11) together with Lemma 6.1 yields 

^ = ^^'^^^ = ^(^,f^) = [J,^] , 
where we recall that a := «(cr). For symmetric spaces, (6.11) is therefore equivalent to 

3 =[a,3], ^ = [J, a] . (6.12) 
In particular ^ = J + a is a constant, as in Remark 6.3 above. 
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6.3 Example: G = 50(3), Q = S'^. 

We work with the conventions defined in Remark 4.5. Choose as anchor point a the North 
pole Bz G S"^ and define the map 

(■)3 : so(3) ^ M, ^ (f2)3 := ri ■ e^. 

The isotropy subgroup SO (3) a = C 5*0(3) corresponds to rotations around e^, and so (3) a 
is identified with M via the map Ae^ i— A. The vector bundle so (3)^ is isomorphic to S*^ x M 
via ^ 

so(3)^ ^ 5^ X M, [A, Ae;] ^ (Ae^, A). (6.13) 

In particular, it follows that the space of reduced variables TS"^ ^Q-5o(3)^ can be identified 
with TS"^ X M. The map a^^ defined in (6.7) becomes 

: TS0{3)/S^ TS^ x M, [A, A] ^ (x, x, {k-^k)^) . (6.14) 

Moreover, the mapping (6.6) is, through the correspondence (6.13), 

52xM^so(3), (x,A)^Ax. (6.15) 

The first-order geodesic equations (6.11) on 5*0(3) become 

J = CTxJ, cr = i X a. 

6.4 Cubics and ballistic curves 

In this section we show that certain types of ballistic curves in a symmetric space Q are 
Riemannian cubics. Recall from Remark 6.2 that a ballistic curve is the projection q{t) G Q 
of a geodesic g{t) G G. Geodesies in G are given by the Lagrange-Poincare equations (6.11). 
As explained in Remark 6.4 these are equivalent to 

j=[a,j], (?=[J,a], (6.16) 

where J := J(g). Recall that for a curve g{t) G G with projection q{t) G Q we defined 
a = [g,A{g)] and J = J{q). As a consequence, the conserved right-invariant velocity is 
TRg-ig = J + 0", where a := i{a) with i defined in (6.6). Moreover, 3{t) = Adg(t) J(0) and 
a(t)=Ad,(t)a(0). 

Theorem 6.5. The projection q{t) E Q of a geodesic g{t) E G is a Riemannian cubic if and 
only if at time t = 

[a, [a, [a, J]]] + [J, [J, [J, a]]] =0. (6.17) 

Proof. Since Ad is a Lie automorphism it follows from 3{t) = Adg(t)J(0) and a{t) = 
Adg(t) (7(0) that (6.17) holds at t = if and only if it holds at all times. Furthermore 
we obtain from (6.16) that 
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Therefore, using also (4.22) and (4.20), 



( J + 2 



0", 



J, J 

0", J 



cr — J, J 
J, J 



a,H„ J 



J,V, (J 



[a, [a, + [j, [j,cr]]]. 



The theorem now follows from equation (4.28) for cubics in symmetric spaces. ■ 

It is clear that (6.17) is satisfied if a = at time t = 0. This leads to geodesies q{t) since 
the projection of a horizontal geodesic is a geodesic. Another class of solutions is given by 
the following corollaries 

Corollary 6.6. Let q{t) E Q he the projection of a geodesic g{t) G G that satisfies [J, a] = 
at time t = 0. Then, q{t) is a Riemannian cubic. 

Corollary 6.7. Let q[t) E Q he the projection of a geodesic g{t) G G that satisfies, at time 
t = 0, 

[3,[3,a]]=ca, [a, [a,3]] = cJ , (6.18) 
for c G M. Then, q{t) is a Riemannian cubic. 

Remark 6.8. For the example G = 5*0(3), Q = S"^, the solutions to equation (6.17) are 
illustrated in Figure 6.1. These solutions are fully described by Corollaries 6.6 and 6.7. 
Namely, by means of the identity 

a X (b X c) = b(a ■ c) — c(a ■ b) 

for vectors in M^, (6.17) is seen to be equivalent to 



J X = 0. 



1 1 2 11-^ m 

Possible solutions are given by J x = 0, or by ||ct|| = J at t = 0. The first case is 
equivalent toJ = OorCT = Oatt = 0, and therefore at all times. This corresponds to trivial 
projected curves x(t) = x(0) G S*^, and to projections of horizontal geodesies, respectively. 

II — ii2 ll^l|2 

We proceed to analyze the alternative solution given by ||cr|| = ||J|| . For a given initial 
velocity x(0) = v G Tx(o)>S'^, the projection x(t) G S"^ describes a constant speed rotation of 
x(0) around the axis 

r2 = J + CT = xxx± ||x|| X. 

Hence, the (constant) rotation axis lies in the plane spanned by J and x, enclosing a 45° 
angle with J. The curve x(t) moves with constant speed along a circle of radius ^ = 



I tI|2 
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Fig. 6.1: Ballistic curves and cubics on the 
sphere. For given initial position and velocity, 
there are two types of trajectories that are, si- 
multaneously, projections of geodesies on the 
rotation group (ballistic curves) and Rieman- 
nian cubics. The two trajectories are shown for 
initial position (1,0,0) and initial velocity par- 
allel to the y-axis. In black a unit speed tra- 
jectory along the equator corresponding to the 
projection of a horizontal geodesic on the ro- 
tation group. The blue curves are the circular 
unit-speed trajectories of radius described 
in Corollary 6.7 and Remark 6.8. 



6.5 Second-order Lagrange— Poincare reduction 

To discuss second-order reduction, we introduce the quotient T^^^G/G^ of T^'^^G by the 
natural action of Ga, whose elements will be denoted by [g,g, 'g] G T'^'^'^G/Ga- We will make 
use of the bundle diffeomorphism a^J^ : T^'^'^G/Ga T^'^^Q Xq 2Qa , 

[g, g, g] ^ (g, q, q) x [g, A{g)] © — [g, A{g)] = (g, q, q) x [g, A{g)] © [g, dtA{g)] , (6.19) 
see [GBHRll]. Here, we defined q{t) := Il{g{t)), where g{t) is a curve representing [g,g, 'g], 

(2) 

that is, {q,q,q) = Tg Il{g,g,g). Denoting the right-invariant velocity by ^ = TRg-ig and 
using the definition (6.2) of A together with (6.1), this becomes 



[9,9,9] ^{q,q,q) x [g,Mg-i\,{i)]®[g,Mg-i\,{i)]. 



(6.20) 



Remark 6.9. In the example G = 5*0(3), Q = S'^, working with the conventions laid 
out in Remark 4.5, the space of reduced variables T^'^^S'^ Xq 2so(3)^ can be identified with 



The map is 



a^J^ : r(2)SO(3)/Si ^ T(2)S^ 



^xR\ 



[A, A, A] ^ (x, X, ±, (A-iA)3, dtiA-'Ah) . (6.21) 



The reduced Lagrangian. For L : T^^^'G — t- M a Ga-invariant Lagrangian we define the 

(2) 

i o a\ . Consider the Lagrangian for Riemannian cubics 



reduced Lagrangian i hy L 

19 



^ ll^^^ll" ^^'-^ yiote the following equalities. 



L{9,g,g) 



D 

(e)Q(^) 

D 



1 



2 

+ 2 

q I 



5 2 



Dt 



^-v,(v,(e))Q 



(6.22) 
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where we used the right-invariance of L, the definition of the normal metric, and part (iii) of 
Proposition 3.1. It follows from (6.20) and (6.22) that the reduced Lagrangian i : T^'^^Q Xq 
2ga — ^ reads 



e{q, q, q, a, a) 



D 
Dt 



(6.23) 



where we recall a := ^cr- The reduced Lagrangian therefore measures the deviations from 
the geodesic Lagrange-Poincare equations (6.11). 

Remark 6.10. In the example G = 5*0(3), Q = S'^, the reduced Lagrangian i : T^'^^S'^ x 

10)2 



^ M IS 



1 2 

(x, X, X, a,&) ^ - II Ax - ax X x||^ + -a 



1.2 1 



I Axjl^ - aDt± ■ (x X x) + Ja^ ||x||J + Ja^ 



where the norm 



is evaluated as the standard Euclidean norm. 



Coupling. The reduced Lagrangian couples the horizontal and vertical parts of the motion 
through the term VgCXg. This explains the absence of a general horizontal lifting property 
for Riemannian cubics studied in Section 5.2. Indeed, let us instead define the Lagrangian 



Lkk : T^^)^ 



with reduced Lagrangian 



{9,9,9) ^ 



D 
Dt 



TU{g) 



(6.24) 



T^^^Q Xq2q,^R, {q,q,q,a,&)^ 



D 
Dt' 



0" 



(6.25) 



The Lagrangian Lj^k belongs to a class of Lagrangians that were studied in [GBHRl 1] as nat- 
ural second-order generalizations of the Kaluza-Klein Lagrangian. The reduced Lagrangians 
£ in (6.23) and Ikk in (6.25) differ by the coupling term V^aQ. The decoupled form of Ikk 
leads to a general horizontal lifting theorem. Namely, any horizontal lift g{t) to G of a cubic 
spline q{t) on Q is a critical point of the action J Lkk dt. 



Lagrange-Poincare equations. We now compute the Lagrange-Poincare equations. Tak- 
ing e-variations and defining Vq := -^g — Vgcrg, we have, for the first term of (6.23), 



D 
Dt 



dt 



,^DD. D ^ 



dt. 



We then compute 
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and 



Dt 

Lemma 6.1 shows that 



^ ' Vsq^Q + (5^)q {q))+R (Sq, q) agiq) - i5dta)Q (q) - Vsg idta)^ 



since 5a := o"e ^ sj'- So we have 

^ ' D \ , f f D 



and 



{5-a)Q (q) = iti6a))Q (g) + a))^ (g) = (F(<5g, a))^ (g), 

7Q (q) jdt = -J^^Q l^jyV,, {F{6q, a))^ (q) J dt 

\Q{Vg,{5dta)Q{q))dt= [\{j{V,),dtSa)=- [' j {dt3{V,),6a) 

Jo Jo 

^{dt3{V,),t{5a) + FiSq,a))=- [ ^ {t;dt3iV,),5a) + {Fjdt3{V,),5q) 

Jo 



where : — )■ TgQ is the transpose of the map Fg-^ : TgQ — )■ q, Ffj^{vq) := F{vq,aq), and 
: — 7- (0a)q IS the transpose of the map ig : {Qa)q (the restriction of (6.6) to the fiber 
{9a)q of Qa at q). 

For the second term of (6.23), we have 



1,1 .,,2 , f. D^. , r 

olkllg„«t = / 1 {(T, —(T ] dt. 



De 



Using the variations 



8(J = —7] - [a, rj] + B{5q, q) E Qa, = -j^^o- - [B{q, 6q), a] G Qa 
and the formula 

|j(«(t))-j(^a(t)) =^^(a(t),g(t)), 
where {^F^ {aq,Vq),rj^ := {^a^, V^^r^g), (see [GBHRIO]) we get the equations 

+ VaJ ■ ^Vq + V{dta)l ■ Vq + R{Vq, q - aQ{q))q 



+ + adt d + {"qdiiiVq), iqB ) = FJ (^^ (v;\ g) ) ' 



I ^ + adt) (^^a + ^^9tJ(K)^ = 0, 
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where we recall that Vg := — Vgcrg G TQ. 

Using (6.10) these equations can be rewritten as 



+ R{VgcyQ, q - aQ{q))q + (z (a + ad^ a + t"^ dtJ (V,))) ^ + FJ (^^ {V^, q)f (6.26) 
— + adtj (a + zj9,J(y,))=0. (6.27) 

These equations are the second-order analogue of (6.11). The left hand side of the first one 
is the equation for Riemannian cubics on Q. Hence the right hand side of (6.26) is the 
obstruction for the projected curve to be a cubic. For symmetric spaces, solutions to (6.26)- 
(6.27) with vanishing obstruction include in particular the curves characterized in Theorem 
5.1, but also the special geodesies on the group that project to the ballistic curves of Section 
6.4. 

Remark 6.11. For G = 5*0(3), Q = S^, the Lagrange-Poincare equations are computed as 



^^gX -I- R yj^^i xj X = cr X + 2(T(Tx + (t(x X x) + 3cr(x X x) 
+ (t(2[x X X + (x X x)"*"] + ||x|p(x X x)) — a(x x x) 



X ■ (x X x) + a||x|p — a = a , 



for a constant a G M. Here we denoted by V"*- = v — x(x ■ v) the orthogonal projection of v 
onto the tangent plane to S"^ at x. 



7 Summary and Outlook 

This paper has investigated Riemannian cubics in object manifolds with normal metrics and 
in particular their relation to Riemannian cubics on the Lie group of transformations. 

Our starting point was the definition, in Section 2, of necessary concepts and a treatment 
of covariant derivatives for normal metrics in Section 3. The derivation of the Euler-Lagrange 
equation for cubics from the viewpoint of normal metrics followed in Section 4. The examples 
of Lie groups with invariant metrics and symmetric spaces were discussed in detail, and the 
relation with equations previously present in the literature was clarified. The new form of 
the equation was seen to lend itself to the analysis of horizontal lifts of cubics, due to the 
appearance of the horizontal generator of curves. 

Section 5 proceeded with this line of investigation by deriving several results about hor- 
izontal lifting properties of cubics. For symmetric spaces a complete characterization was 
achieved of the cubics that can be lifted horizontally to cubics on the group of isometries. In 
rank-one symmetric spaces this selects geodesies composed with cubic polynomials in time. 
The section continued with a treatment of the corresponding question in the context of Rie- 
mannian submersions. In Section 6 certain non-horizontal geodesies on the group were shown 
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to project to cubics in the object manifold. A complete characterization of such geodesies 
was given in the sense of Theorem 6.5. For the unit sphere acted on by the rotation group 
the corresponding projections were seen to be the circles of radius 1/a/2. A discussion of 
Lagrange-Poincare reduction of cubics led to reduced equations that identified the obstruc- 
tion for projections of cubics to be cubics in the object manifold. 

Future research should seek to extend the results found here for horizontal lifts of cubics 
in symmetric spaces to more general situations. A first step in this direction has been 
taken in the context of Riemannian submersions in Section 5.2. However, a more complete 
characterization of the cubics that lift horizontally to cubics may provide a wider class of 
applications. A similar remark holds for the analysis of the special ballistic curves that satisfy 
the equation for cubics. All of these problems are tightly linked to the obstruction term in 
the Lagrange-Poincare equations, particularly on the right hand side of (6.26). Hence, one of 
the main tasks ahead is to deepen the understanding of that obstruction term, and thereby 
determine additional situations in which it vanishes. 
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